Path Integral Ground State study of 2D solid "^He 
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We have studied a two-dimensional triangular commensurate crystal of *He with the exact T = K 
Path Integral Ground State (PIGS) Monte Carlo method. We have projected onto the true ground 
state both a Jastrow-Nosanow wave function, in which equilibrium positions are explicitly given 
and no Bose-Einstein (BEC) is present, and a translationally invariant shadow wave function, in 
which the solid phase emerges through a spontaneously broken symmetry process and it has BEC. 
We find a remarkable convergence to the same properties, both the diagonal ones as well as the 
off-diagonal one-body density matrix pi. This supplies a strong evidence that no variational bias 
are present in the PIGS method. We find no BEC in the commensurate 2D *He crystal at T = K, 
p\ shows an exponential decay in the large distance range. The structure found in p\ is due to 
virtual vacancy-interstitial pairs and this shows up in the momentum distribution. 

PACS numbers: 67.80.-s 
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Quantum Monte Carlo methods have provided a very 
powerful tool in exploring the physics of strongly inter- 
acting many-body quantum systems. As far as the prop- 
erties of Bose fluids and solids are concerned. Path Inte- 
gral Monte Carlo (PIMC) methods have been proved to 
evaluate "exact" and "unbiased" expectation values on 
the thermal equilibrium state at finite temperature 1']; 
"exact" means that the obtained results are the true ex- 
pectation values within the statistical error, while "unbi- 
ased" means that the only required input is the inter- 
atomic potential. At zero temperature, different "ex- 
act" techniques are available, as the Green's function 
Monte Carlo (GFMC), the Diffusion Monte Carlo [3] 
(DMC), the Reptation Monte Carlo [i] or the Path Inte- 
gral Ground State (PIGS) methods. All such meth- 
ods, though "exact" in suitable limits, rely on variational 
models for the ground state wave function (wf) of the 
system. These variational models play a relevant role 
in the importance sampling employed by these methods; 
although, in principle, the final results should not be af- 
fected by the particular choice of the trial wf, in practice, 
the possibility that some bias could survive, especially 
in systems with complex broken symmetries like in the 
solid phase, has not been established. In this paper we 
will show the robustness of the results given by the PIGS 
method with respect to the choice of the variational wf 
by studying a model for *He in two dimensions in the 
solid phase. In fact, we find the same physical properties 
of the system, within the statistical noise, starting from 
two radically different wfs. One, a Jastrow-Nosanow Q 
wf (JNWF), has built in the crystal lattice via Gaus- 
sian localization factors, it is not Bose-symmetric and 
has no Bose-Einstein condensation (BEC). The other, a 
shadow wf (SWF), is translationally invariant, with 
crystalline order arising as spontaneous broken symme- 
try, and it has BEC [l^ . We have chosen to study 
this model in 2D for different reasons. The reduced 
dimensionality allows us to study correlations to much 
larger distances than in 3D, even more than 20 lattice 



parameters. Fluctuations are expected to be stronger in 
2D so this is a more stringent test for convergence given 
the different symmetry properties of the starting varia- 
tional wfs. Finally, the 2D system is a model which is 
relevant for adsorbed "^He atoms on a smooth substrate 
like graphite. 

With the present study we address also the im por tant 
question of the supersolid state of matter 11, l^, 13, 14, 



151 . Il6j . PIMC computations give strong evidence that in 
a 3D commensurate solid ^He (number of atoms equal to 
the number of lattice sites) there is no superfluid response 



and no BEC [17|, [1J| . The PIMC computations are at 



finite temperature and the lowest T is of order 0.1 K 
which is above the experimental transition temperature 
of crystals of good quality 12]. By computing the density 
matrix pi{r,r') of crystalline *He at T = K, we find 
that there is no BEC in a 2D crystal. Preliminary results 
indicate that this is true also in 3D. 

Dealing with low temperature properties, "^He atoms 
are described as structureless zero-spin bosons, interact- 
ing through a realistic two-body potential, that we as- 
sume to be the IIFDHE2 Aziz potential [T?| . The aim of 
the PIGS method is to improve a variationally optimized 
trial wf by constructing a path in the Hilbert space of the 
system which connects the given wf to the true ground 
state of the system; during this "path" , the correct corre- 
lations among the particles arise through the "imaginary 
time evolution operator" e~'^^ , where H is the Hamilto- 
nian operator. Being (jj a trial wf with non-zero overlap 
with the exact ground state "001 this V'o is obtained as the 
T — > oo limit of tpT = e~'^^ (j) suitably normalized. This 
tpr can be written analytically by discretizing the path in 
imaginary time and exploiting the factorization property 
q-(ti+t2)h _ g-TiJ?g-T2ff ^ this way, i/'r turns out to 
be expressed in term of convolution integrals which in- 
volve the "imaginary time propagator" {R\e~^'^^\R') for 
a St, that can be small enough such that very accurate 
approximants are known [J [2o| . This maps the quan- 
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turn system into a classical system of open polymers Q . 
An appealing feature peculiar to the PIGS method is 
that, in ip^, the variational ansatz acts only as a starting 
point, while the full path in imaginary time is governed 
by e""^^, which depends only on the Hamiltonian oper- 
ator. 

As a trial wf, we used both a JNWF and a SWF. The 
JNWF is written as the product of two-body correla- 
tions and of Gaussian one-body terms which localize the 
particles around the assumed lattice positions j6j. In 
the SWF, beyond the explicit two-body factors, addi- 
tional correlations are introduced via auxiliary (shadow) 
variables which arc integrated out 0. Nowadays, SWF 
gives the best variational description of solid and liquid 
^He fill. In both the cases, the variational parameters 
have been chosen to minimize the expectation value of 
the Hamiltonian operator. In what follows, we will refer 
to PIGS when we deal with the projection of the JNWF, 
and to SPIGS ^ (Shadow Path Integral Ground State) 
when we project the SWF. 

Because of the Bose statistics obeyed by the atoms, 
when using ili^ as an approximation of the true ground 
state TpQ, one has, in principle, also to account for permu- 
tations in the propagator {R\e^^'^^\R') 1, 5]. Permuta- 
tion moves are necessary when the JNWF, which is not 
Bose-symmetric, is used. On the other hand permutation 
moves are not necessary whenever the trial wf in ?/;,- is 
already Bose-symmetric, as the SWF. However, also for 
SPIGS, adding permutation moves turns out to be useful 
in improving the efficiency and the ergodicity of the sam- 
pling, mainly in reaching the large-distance range of pi. 
In our algorithm we have introduced two different permu- 
tation samplings: cycles of inter-particles exchanges and 
swap moves. The first, which may involve an arbitrary 
number of particles, are described in detail in Ref. l24l . In 
a PIGS or SPIGS calculation of pi(r, r'), the efficiency 
can be further improved by introducing particular two- 
particles permutations cycles involving always one of the 
two positions r and r': the swap moves fl^. These moves 
improve very much the efficiency of the computation of 
pi and their acceptation rate is remarkably high: in the 
present 2D system, by using a staging 25] method to sam- 
ple the free (kinetic) part of the imaginary time propa- 
gator, we have found an acceptation rate for this swap 
move which is nearly 15% in the liquid phase and in the 
solid phase at densities close to the melting. 

The 2D "^He system phase diagram is known from ac- 
curate finite temperature PIMC simulatio ns [26l| ; at zero 
temperature both DMC [27|] and GFMC (isfhave been 
used to investigate its properties mainly for the liquid 
phase. We have performed SPIGS and PIGS simula- 
tions of a 2D ■^He commensurate triangular crystal at 
p — 0.0765 A^^, slightly above the melting density. In or- 
der to control the reliability of our results we have tested 
their dependence on both the "projection time" t and the 
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Figure 1: Energy per particle E/N and Debye- Waller fac- 
tor e~^^ as functions of the projection time r obtained from 
PIGS and SPIGS for a commensurate 2D '*He crystal with 
TV — 120 particles a,t p — 0.0765 A~^. Error bars are smaller 
than the used symbols. The dotted lines indicate the conver- 
gence values eo = 1-308 ± 0.002 K and 6"^^ = 0.214 ± 0.002, 
solid lines are exponential fits to guide the eye. 



"time step" St. For a fixed value of t = 0.075 K~^, we 
have done calculations with St = 1/40, 1/80, 1/160 and 
1/320 and we have used the pair-product approxi- 
mation [l[ for the imaginary time propagator. Reducing 
St below 1/40 affects only marginally the results; 
for example, by using 1/320 K~^, the obtained energy is 
only 1% lower then the one with 1/40 K^^. So we have 
adopted in most of the computations St = 1/40 as 
a reasonable compromise between accuracy and compu- 
tational effort. These tests provide also a robust check 
of the ergodicity of the sampling algorithm since a lower 
value of St for given r means a bigger number of small 
time projections (convolution integrals) so that one is 
dealing with polymers of increasing length. We have 
then increased t till convergence in the results has been 
achieved. 

Diagonal properties, like the energy, have been com- 
puted in a box which hosts N = 120 particles with peri- 
odic boundary conditions. In Fig. [1] we give the energy 
per particle as a function of t both for SPIGS and PIGS: 
aX T — 0.075 K"-'^ the SPIGS result is already converged 
to the value E/N = 1.309±0.002 K, while with the PIGS 
one reaches convergence at r = 1.575 K"-'^, where the 
energy takes the value E/N = 1.311 ± 0.002 K. The po- 
tential energy values are respectively —11.053 ± 0.004 K 
and -11.048 ± 0.006 K. From Fig. [H it is evident that 
the true ground state value is reached much more quickly 
(i.e. with a lower number of small-time projections) with 
SPIGS. The overlap between the trial wf and the true 
ground state plays a crucial role in the convergence: a 
large overlap ensures that i/'r is a good approximation 
even for not too large r. It has been shown [2§| that 
|(?/'o|?!')P = (e~''')^ where 7 corresponds to the area of 
the region between E{t)/N and its convergence value eq. 
From our results, e~'^ turns out to be 99.8% for the SWF 
and 97.9% for the JNWF. From the peaks in the static 
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Figure 2: One body density matrix pi{f—f") computed along 
the nearest-neighbor direction in a 2D commensurate ''He 
crystal at the density p = 0.0765 for different r values. 
Notice the completely different behavior for the two initial 
wfs and the convergence to the same pi. 



structure factor we have extracted the Debye- Waller fac- 
tor; we have found that both SPIGS and PIGS converge 
to the same value e''^^ = 0.214 ± 0.002 (Fig. [1]) and, 
as in the energy case, SPIGS shows a faster convergence. 
It is important to note that the energy convergence does 
not allow to deduce a priori the convergence of other 
physical properties: the convergence must be checked in- 
dependently for each observable. 

In order to study whether the 2D commensurate solid 
''He has BEG, we have computed the one-body density 
matrix pi(r, r '). A non-zero limit of pi for |r — f'\ ~^ oo 
implies BEG and the Fourier transformation of pi gives 
the momentum distribution n^. We have computed pi 
in a 2D commensurate crystal at p = 0.0765 in a 
simulation box with = 240 particles (this allows us 
to explore distances up to about 39 A). In our calcula- 
tions we have sampled pi(r, r') in two different ways. 
In the first one we force r — r' to lie on a nearest- 
neighbor direction, while, in the second r and r ' are al- 
lowed to explore the full plane; we have found perfectly 
consistent results. In Fig. [5] we show the convergence of 
pi computed along the nearest-neighbor direction with 
SPIGS and PIGS for increasing projection time. The pi 
given by SWF has a non-zero limit at large distance, 
i.e. there is BEG as in 3D When projecting from 
SWF the large distance plateau decreases for increasing 
T until, at T = 0.775 K"-', it disappears up to the dis- 
tance allowed by the simulation box. At convergence, pi 
displays an exponential decay, with a correlation length 
A ~ 2.75 A, with superimposed a small modulation re- 
flecting the crystal symmetry. The pi given by the JNWF 
is essentially a Gaussian. Under projection, exchanges 
among the atoms greatly extended the range of pi . It is 
remarkable that at t = 0.775 K"'^ the pi obtained start- 
ing from such different wfs coincide within the statistical 
uncertainty. Similar conclusions are reached also from 
the computation of pi in the whole plane. In Fig. [3^ we 



show pi obtained both with SPIGS {x — x' > Q range) 
and PIGS {x — x' < range) with the projection time 
T = 0.775 K^'^. Again, the results are indistinguishable 
within the statistical error. We have looked also for fi- 
nite size effects by considering system at the same density 
but with different particles number [N — 180, 240 and 
480) and we have found no appreciable differences. With 
PIGS we have studied pi up to a distance of about 60 A, 
finding the same exponential decay. 

We conclude that a 2D commensurate ''He crystal has 
no BEG, or, more precisely, the condensate fraction, if 
any, is below 10^^. We notice that exchange processes 
are rather significant in the system, so that the range of 
pi is significantly larger than the size of the unit cell. This 
manifests itself in the momentum distribution (Fig. [3)3) 
as deviation An^ of from the Gaussian correspond- 
ing to the kinetic energy. There is an excess of particles 
at low momenta up to fc = 1.2 A^'^ whereas there is a 
deficit in the A:-space region around k = 1.6 A~^. The 
positions of the bumps of pi over the exponential decay 
suggest that the extension of pi beyond the unit cell can 
be interpreted as a signal of the appearence of vacancy- 
interstitial pairs (VIP): their position respect to the ori- 
gin (r — r') correspond to interstitial positions whereas 
the distance between the positions of two neighboring 
pumps is equal to the lattice parameter. In the case of 
the SWF these VIP are unbound allowing for a finite 
BEG [§], whereas these pairs are bounded in the exact 
i/'o- This might be due to the presence of some long range 
correlations in -00, possibly caused by the zero-point mo- 
tion of transverse phonons, which are absent in the SWF. 
Since the ground state wf reflects the zero-point motion 
of any excitation in the system [lOl , it is tempting to sug- 
gest that the evidence of VIP in pi is the manifestation 
of a new kind of excitation in the quantum solid beyond 
the phonon excitations. 

Our results that the 2D commensurate cr ysta l has no 
BEG at T = K and the similar results [13, [3 ob- 
tained in 3D at flnite T make plausible that also in 3D 
a commensurate crystal has no BEG at T = K and, 
in fact, we have some preliminary results for pi in 3D 
supporting this notion. Experimentally it is established 
that defects and '^He impurities have an important role on 
the non-classical rotational inertia (NGRI) of ^He crystal 
and dislocations are suspected to have a relevant role [l^ , 
at least in the case of good quality crystals. It is an open 
question what happens in a solid with less and less de- 
fects: will any NGRI effect go away? Our result suggests 
that it should go away unless some disorder is present 
even in the ground state as an intrinsic property of ■00- 
A key question is therefore to establish whether zero- 
point defects are present in the ground state of a quantum 
solid. We believe that the presence of zero-point defects 
in solid ^He is still an open question [2^. Since a SWF, 
which has BEG, is the exact ground state of a suitable 
Hamiltonian, the interesting question is for which class 
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Figure 3: (a) One-body density matrix pi(r — r') computed in a 2D commensurate "^He crystal with 224 particles at the density 
p = 0.0765 A~^, in the j: — a;' > (a; — a;' < 0) range the SPIGS (PIGS) results are displayed, (b) Deviation of the momentum 

distribution Ang from the Gaussian distribution nj: = exp ^— gjfy) ^ ; obtained by Fourier transforming the difference 

between the SPIGS pi and the Gaussian distribution corresponding to the kinetic energy (T). 



of interatomic potentials does the commensurate crystal 
have BEC and supersolidity? 

The convergence in the results obtained starting from 
radically different wfs is a remarkable result because sup- 
plies strong evidence for the absence of any variational 
bias in the PIGS method. Moreover, even if not shown in 
this letter, we have also obtained convergence of diagonal 
properties, such as the static structure factor and the ra- 
dial distribution function, by projecting a simple Jastrow 
wf which has just the minimal information on the short 
range behavior and displays no crystalline order at the 
considered density. 
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